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КОМПОЗИЦИОННЫЙ СИНТЕЗ КВАЗИОПТИМАЛЬНЫХ 
ПО БЫСТРОДЕИСТВИЮ СИСТЕМ УПРАВЛЕНИЯ 
ВЫСОКОГО ПОРЯДКА 


Предложен аппроксимационный метод нахождения квазиоптимальных по быстро- 
действию управлений гладкими стационарными системами, удовлетворяющими 
ряду структурных ограничений. Построена методика синтеза законов реализации 
обратных связей, реализующих квазиоптимальное по быстродействию управление 
нелинейными системами. Особенностью методики является то, что она позволяет 
использовать результаты синтеза системы низкого порядка для построения закона 
управления системой более высокого порядка. Это позволяет реализовать компози- 
ционный синтез САУ теоретически любого порядка на основе эффективного и на- 
дежного результата, полученного для динамической системы первого порядка. 
Ключевые слова: динамическая система, закон управления, быстродействие, 
квазиоптимальность, синтез, аппроксимация, композиционный подход. 


Введение. Проблема оптимизации законов управления является одной из 
важнейших проблем современной теории автоматического управления 
(СТАУ). Однако для задач оптимального управления характерно, что их 
аналитическое решение удается получить лишь в редких случаях. В связи с 
этим в СТАУ разработаны различные способы нахождения аппроксимаци- 
онного решения таких задач. Результаты подобного решения называют за- 
конами квазиоптимального управления. Большинство методов квазиопти- 
мизации быстродействия направлено на аппроксимацию поверхности 
переключения или на численные методы нахождения моментов переключе- 
ния [1]. 

В настоящее время задача синтеза систем, оптимальных по быстро- 
действию, является одной из наиболее актуальных задач теории оптималь- 
ного управления. Время регулирования входит в число основных каче- 
ственных характеристик систем автоматического управления. Для многих 
технических систем уменьшение времени регулирования, т.е. повышение 
их быстродействия, имеет большое практическое значение. В статье [2] 
предложен, а в работах [3-5] развит новый подход к построению аналити- 
ческих моделей ё -квазиоптимального быстродействия для задач управле- 
ния. Его сущность состоит в статической аппроксимации правой части ма- 
тематической модели (ММ) оптимальной системы, а преимущество заклю- 
чается в том, что качество управления, оцениваемое временем регулирова- 
ния, задается в модели параметрически. Варьирование параметров # поз- 
воляет изменять степень квазиоптимальности решения. 

Опираясь на введенное в [2] понятие квазиоптимальности по вы- 
бранному критерию, можно предложить и другие, отличные от сформули- 
рованного в [2, 3], подходы к аналитическому нахождению квазиоптималь- 
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ных по быстродействию управлений. Один из таких методов, ориентиро- 
ванный на приближение быстродействия проектируемой системы к соот- 
ветствующему показателю уже синтезированной квазиоптимальной систе- 
мы, но более низкого порядка, предлагается в данной работе. 

Постановка задачи. Пусть управляемая система И -го порядка описыва- 


ется системой дифференциальных уравнений т.н. треугольного вида, т.е. 
ие не, (1) 


1 


где функции Х. - дифференцируемы по всем переменным х,,х,,...х., а 


у | 
т й 


Требуется найти управление и= и(х,х,,...х,), переводящее 
0 0 0 
объект из начального состояния (х, ‚Х,...х,„) в начало координат за близ- 


кое к минимальному время ор при ограничениях на все производные 
переменной состояния и на управление 
51 п 1; 





И и. (2) 


Композиционно-динамический подход к структурно-параметри- 
ческой оптимизации свойств управляемой системы. Подход строится 
на предположении, что известно квазиоптимальное по быстродействию 


| 





Ех, 


т 


КуорЕ К К К 
управление и„””[х), являющееся гладкой функцией состояний 


Т 
х= [х,х,,...х„) , для системы (п- 1) -го порядка 


*^“п-1 
И ЕЕ (3) 
у которой первые и-2 функции правых частей, начальные условия 


00 0 
Х >» ›- + М1 И ограничения на переменные СОСТОЯНИЯ 





| |= х„з;Ё Ъи 2 совпадают с аналогичными составляющими исходной 


1 





системы (1), а на управление и,„, наложено ограничение 


(4) 


Поскольку повышение порядка динамической системы (ДС) без из- 
менения ее внутренней динамики (взаимосвязей по переменным состоя- 
ния) приводит к повышению ее инерционности, оценку быстродействия ДС 
(3) можно рассматривать как естественное ограничение для быстродей- 
ствия системы (1). Тогда закон движения ДС (3) можно рассматривать как 
динамический эталон поведения системы (1). 


Однако в (1) правая часть (п- 1) -го уравнения формируется по за- 








и, | а 


кону 
ЛЬ». , (5) 
тогда как для квазиоптимального движения ее структура должна иметь вид 
КуорЕ 
и О : (6) 


Трансформация выражения (5) в форму (6) в управляемой ДС мо- 
жет быть осуществлена только управляющим воздействием, а его синтез 
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может быть осуществлен на основе формирования в управляемой ДС инва- 
риантного многообразия 
КуорЕ 


а и = 0 (7) 


т 
которое можно рассматривать как предел Пте х 1 }) = 0 некоторой гипо- 


тетической переменной, представляющей собой ошибку аппроксимации 
свойств эталонной системы, фактическим законом изменения текущих зна- 
чений переменных состояния 
ь. Куор! 

В оне ыы | (8) 

Для обеспечения автоматической системе таких свойств, чтобы при 
квазоптимальном движении ее (и- 1) -х координат ошибка (8) стремилась 
к нулю, при выполнии условия (7) необходимо связать ее движение требо- 
ванием асимптотической устойчивости по макропеременнойе( х,х,,...,х, | 


. Этого можно, например, достичь, потребовав от Е свойства 





е(х) = - 380. (9) 


где & - параметр настройки решения задачи синтеза [6]. 

Наложение дополнительных условий на связи управляемой системы 
возможно только через движение И -й координаты, поскольку остальные 
связаны постулированным свойством гипотетической подсистемы (3). В 


свою очередь переменная х, (1 зависит от управления и( х) в (1). 


Раскрыть механизм реализации указанной связи можно, если про- 
дифференцировать выражение (8) по времени: 
-1 
п-1 п-1 
ой о, ор 9е 9. 
е= ее ое и. 
ри дх, дх, ты дх, дх, 





(10) 


и подставить в уравнение (9) как (8), так и результат дифференцирования 
(10). 
Это позволяет получить искомый закон квазиоптимального управ- 





ления 
-1 
е п-1 де 9 
Е) Е" 1.— . 9-1 ы т: (11) 
сн 9х, 0х, 
или, после преобразований, 
В 1 АуорЕ _ = 
ие: | 6 В М Ли, п р [мк И а :а У, (12) 
ё 21 9х, 0х, 


Результат (12) показывает, что квазиоптимальное по быстродей- 
ствию управление системы И- го порядка можно находить через квазиоп- 
тимальное управление системы (и- 1)- го порядка. Таким образом, по 
принципу математической индукции теоретически можно получить анали- 
тическое решение поставленной задачи для управляемой системы любого 
порядка, опираясь, например, на модель квазиоптимального по быстродей- 
ствию управления первого порядка, предложенную и исследованную в [2- 
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5] как на эталонную. Условие асимптотической устойчивости предложенно- 
го метода также можно доказать математической индукцией, так как для 
системы (п-1)-го порядка (3) существует квазиоптимальный закон быстро- 
действия, обеспечивающий глобальную асимптотическую устойчивость. 
Следовательно, существует — соответствующая функция Ляпунова 


Г(х ,х,,...х„1). Для системы п-го порядка в качестве функции Ляпунова 
можно выбрать функцию 


1 
о О ее (аьжь,..-д,). (13) 


Тогда из (9), (13) можно получить 
Инне 0’ 


(14) 


1 
(хх, х,,...х,) = И(х,х,,...Х.1)- 2-е (х,х,,...х,)< 0\х 0. 
ё 


Формула (14) показывает, что функция И’(х,,х.,...х,) действительно яв- 


ляется функцией Ляпунова для исходной системы (1) — асимптотическая 
устойчивость системы квазиоптимального быстродействия доказана. 

Однако сложность аналитических преобразований при существен- 
ной нелинейности как системных, так и аппроксимирующих функций, вхо- 
дящих в многообразие (7), ограничивает фактический максимальный поря- 
док синтезируемой системы. В многомерном случае, когда управление яв- 
ляется векторным, но имеется возможность декомпозировать исходную си- 
стему на треугольные подсистемы, предложенную методику можно приме- 
нить для каждой подсистемы в отдельности. 

Существо и алгоритм реализации разработанного и описанного в 
данной работе метода наглядно иллюстрирует приводимый ниже де- 
монстрационный пример. 

Иллюстрационный пример. Синтез квазиоптимального по быстро- 
действию управления движением спутника в центральном поле 
сил на круговой орбите. Дифференциальные уравнения динамики плос- 
кого движения спутника в поле земного тяготения, согласно [7], имеют вид 


: у : У У 1 

ЖЕ + защ, = аа = 225+ —и,, (15) 
У | У 

где У =У, М; у, - постоянная все- 


мирного тяготения, Л - масса Земли; 
у = А- Расстояние от ее центра до 


спутника; У: = © - угол поворота при 


движении спутника в полярных коор- 
динатах относительно земли (рис.1), 
у, -скорость его движения вдоль ра- 


диус-вектора А; У.- угловая ско- 
рость вращения спутника. 

Требуется найти управления и, и., 
переводящие спутник из некоторого 





Рис. 1. Схема движения 
спутника вокруг земли 
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0 0 0 0 
начального состояния и. на заданную круговую орбиту 


У = К, у, = 0, у; = уаг, у =@ "| за время, близкое к минимально воз- 


можному. 
Из результатов, полученных в работе [21, известно, что квазиопти- 
мальное по быстродействию управление с ограничением и" для системы 
КуорЕ „ ^ ПИ х т 
первого порядка имеет вид и) ыы РВ и” = сои8Ё> 0. Поскольку 
Ухо +ё 


угол вращения у, = @ может принимать произвольные значения, то можно 


применить изложенную методику для двух подсистем: Я = У>, 


: у Е у. 1 
У =- я" ху: ии У. = р. ие в отдельности. Это иллюстри- 
1 1 1 


руется приводимыми ниже выкладками. 
Для первой подсистемы 














и! в м Ол Й К) : й бе в. 
1 до НИ 5 ЧЕЛА 
[»-2,) +Е, Е 
откуда 
и= у+ и (- №) 1, у у уе" 
1 2 5 2 1.74 р 1,5 
[у- №) +27 И. (»- А) +21 
Для второй подсистемы 
КуорЕ _ - и? (у.-б.) 
й = 2 2! (16) 
[4-й +2, 
откуда 
и» -б 
и, = ЕВ 2... (17) 


2 2 
(5. -6 и +8, 
Нетрудно заметить, что квазиоптимальные аппроксимации управле- 
ний обеспечат также ограничения на производные |5; | < м”, |у4|< и". 


Результаты имитационного моделирования синтезированной САУ 
на ПЭВМ. На рис.2 изображены переходные процессы для моделирования 
координат управляемого законами (17) и (18) движения спутника со следу- 


ющими параметрами: у = 3, 0 = 110 = 5,0 =3; №=2; б@0=4; 
и" = м = 1; Е = 0.Цу = 100. На рис. 3 приведены помеченные индекса- 
ми траектории движения спутника из различных начальных состояний: 
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0 0 0 И 0 0 0 "АТ 
Й = № = 4, у = 1, уз = 5,4 = 3}; Г. = № = 3, у2 = Ъ уз = 0.5, уд = 4}; 


0 0 | 
Г: = [= 2.5, 0 = 0.5, уз = 4, уд - 3}; 
0 0 
== 4.79.59 = 2,0 = 2]. 












































Рис.2. Переходные процессы Рис.3. Траектории движения спутника 
квазиоптимального движения спутника из разных начальных состояний 


Выводы. На основе изложенного в статье метода можно аналитически 
синтезировать квазиоптимальные по быстродействию законы управления 
для нелинейных объектов любого порядка при ограничениях как на фазо- 
вые координаты, так и на само управление. При этом результат решения 
задачи синтеза может быть представлен аналитической функцией перемен- 
ных состояния, процесс нахождения которого можно реализовать, в 
частности, и в виде программного алгоритма для ЭВМ (например, с помо- 
щью пакета 5утБойс Ма ТооБох в МаМаБ). 
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